Superfluidity of an interacting trapped quasi-2D Bose gas 
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We investigate the harmonically trapped interacting Bose gas in a quasi-2D geometry using the 
classical field method. The system exhibits quasi-long-range order and non-classical rotational 
inertia at temperatures below the Berezinskii-Kosterlitz-Thouless cross-over to the superfluid state. 
In particular, we compute the scissors-mode oscillation frequencies and find that the irrotational 
mode changes its frequency as the temperature is sweeped across the cross-over thus providing 
microscopic evidence for the emergence of superfluidity. 
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I. INTRODUCTION 

Superfluids are peculiar states of matter in which, at 
the cost of losing part of their individuality, particles 
gain the ability of cooperative lossless motion The 
occurrence of superfluidity can be attributed to the for- 
mation of certain non-local correlations within the sys- 
tem. For two-dimensional (2D) systems, Mermin and 
Wagner showed that Heisenberg models can be neither 
ferromagnetic or anti-ferromagnetic at finite tempera- 
tures [2j. Hohenberg further ruled out the existence of 
long-range ordering in 2D Bose and Fermi systems 
These rigorous results imply that Bose-Einstein conden- 
sation (BEC) does not exist at any finite temperature in 
uniform, interacting 2D systems in the thermodynamic 
limit, since spontaneous long-range ordering is prevented 
by long-wavelength fluctuations. Hence one might ex- 
pect that conventional superfluidity would not occur in 
2D systems. 

Nevertheless, a different path to superfluid behavior is 
possible in 2D systems. At low temperatures quasi-long- 
range correlations may form with an associated power- 
law decay that eventually reaches zero instead of extend- 
ing throughout the system. It was theoretically shown by 
Berezinskii |4( and by Kosterlitz and Thouless 5} (BKT) 
that in 2D a transition to a superfluid state may occur 
at finite temperatures. Qualitatively the physics of such 
low-tempcrature 2D system is conveniently described us- 
ing the notion of vortex-antivortex pairs (VAPs). At 
temperatures below the BKT transition, long-wavelength 
fluctuations destroy true long-range order and yield spon- 
taneous creation and annihilation of bound VAPs at the 
boundaries of the local domains of the resulting "quasi- 
condensate" Q . These locally coherent blocks contribute 
to the power-law decay of the two-point correlation func- 
tion resulting in a superfluid response of the system. 
On increasing the temperature, fluctuations increase and 
VAPs unbind at the BKT transition. The breaking of 



VAPs results in the proliferation of free vortices and an 
exponentially decaying correlation length, and hence the 
system loses its superfluid properties. This was further 
quantified by Nelson and Kosterlitz who predicted a uni- 
versal jump in the superfluid density at the critical point, 
which may be used to empirically detect the BKT tran- 
sition 0. 

In a trapped ultra-cold Bose gas the situation is rather 
complicated due to the inhomogeneity arising from the 
confining potential. Bagnato and Kleppner showed that a 
trappedj ideal Bose gas undergoes BEC at finite temper- 
atures [8j , implying that in principle a trapped, interact- 
ing Bose gas may exist in a coherent BEC phase and/or 
in a BKT-type phase. The low temperature structure 
of a real 2D Bose gas in a trap has therefore attracted 
a fair amount of theoretical discussion in the recent lit- 
erature and it has been debated whether the superfluid 
transition in such systems is of BEC-type or BKT-type 

Experimentally, the superfluid BKT transition in a 
bulk was first realized in liquid helium thin films by 
Bishop and Reppy [2l|. Resnick et al. reported an ob- 
servation of the transition in superconducting Joseph- 
son junction arrays [22j], followed by Safonov et al. who 
measured a kink in the three-body loss rate in spin- 
polarized hydrogen [23} . These experiments relied on in- 
direct methods of observation where as in the trapped 
atomic gases VAPs and their dynamics can be directly 
imaged. The quasi-2D regime in trapped quantum degen- 
erate gases has been approached experimentally by using 
tight axial confinement with the aid of optical potentials 
and through centrif uga l exp ansion in rapidly rotated con- 
densates [24|, |25|, |26|, |27J, |28| . However, probing the details 
of the quasicondensation transition has onl y re cently be- 
come experimentally accessible (29l. l30l. l3ll. |32| . 

Indeed, Dalibard and co-workers observed phase de- 
fects in the interference patterns of multiple quasi-2D 
gases trapped in the valleys of an optical lattice [28| . 
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These phase defects were evidently caused by unbound 
free vortices [l6|, [Hj]. Further observations on spatial 
phase correlations in the system provided evidence for the 
cross-over between the BKT quasicondensate and nor- 
mal state [29J, l2Q| • While the observed correlations were 
shown to be consistent with the BKT quasicondensation 
picture, the question whether the system is superfluid or 
not still remained. Recently an observation of the BKT 
cross-over has been achieved in a 2D lattice of Josephson- 
coupled BECs [3l[ and in a single 2D dipole trap [HJ . 

The interplay between interactions and inhomogeneous 
effects arising from the trapping potential have made 
theoretical predictions for the low temperature phases 
of dilute, atomic Bose gas difficult. The main point we 
address in this paper is the question of superfluidity in 
this system at low temperatures. There is not a single 
observable that categorically defines superfluidity and we 
present results of microscopic calculations for a variety of 
observables, including off-diagonal long-range order (con- 
densate), fluctuations, scissor mode dynamics and pres- 
ence of vortices. From the combined analysis of these 
quantities we are able to infer a cross-over temperature, 
T s f, below which the system exhibits superfluidity. 

The paper is organized as follows: In Sec. II we discuss 
various measurables useful in examining the superfluid 
properties of our system. Our computational approach 
is explained in Sec. Ill, and the results are presented in 
Sec. IV followed by the concluding remarks in Sec. V. 

II. EVIDENCE OF SUPERFLUIDITY 

A substance which has the ability to flow without dis- 
sipation is superfluid. Although the difference between a 
superfluid and a classical fluid may seem intuitively clear, 
it is difficult, if not impossible, to find a single universal 
definition for superfluidity against which any material 
could be tested. Indeed, the complete description of su- 
perfluidity is not a single feature but a complex of phe- 
nomena [33j |. In an interacting 3D atomic gas, the forma- 
tion of BEC is essentially equivalent to the emergence of 
a macroscopic wave function, which inherently exhibits 
long-range order throughout the system. Furthermore, 
the system attains finite superfluid fraction at the BEC 
transition. The situation is more subtle in 2D where 
the condensation process is plagued by long-wavelength 
phase fluctuations. In the following, we introduce mea- 
surables relevant for providing evidence of superfluidity 
in a quasi-two-dimensional sample of trapped ultra-cold 
atoms and apply these definitions to discuss the super- 
fluidity of quasicondensates. 

A. Role of quantized vortices 

Quantized vortices are the hallmark of superfluids. 
The flow v s (r, t) = H/m\7ip(r, t) of a superfluid described 
by a macroscopic wave function with a phase <^>(r, t) obeys 



the condition of irrotationality, V x v s (r, t) — 0, and the 
Onsager-Feynman quantization of circulation, 

<j)v s (r,t) • dl = K2nh/m, (1) 

where n is an integer and m is the mass of an atom. 
Therefore rotation (k ^ 0) is only possible around the 
phase singularity at the core of a quantized vortex where 
the superfluid density vanishes. If normal fluid is present, 
it occupies the volume in the vortex core. The quantized 
vortices can therefore be seen to play a two-fold role in 
superfluid systems. On the one hand their coherent role 
is vital in enabling superfluids to rotate, while on the 
other hand they can be viewed as (topological) defects 
that are a source of incoherence causing a reduction in 
the superfluid fraction. 

It is worth noting that the observation of a vortex alone 
is not a sufficient criterion from which superfluidity may 
be deduced. For instance, while a persistent vortex in a 
zero temperature condensate is readily accepted as proof 
of superfluidity, in a 2D system near BKT cross-over an 
observation of a transient isolated vortex would be more 
likely to signify unbinding of VAPs and hence loss of 
superfluidity. Therefore, special care must be taken in 
the treatment of a situation where both spontaneously 
(thermally activated) and actively (by external rotation) 
created vortices may exist simultaneously. A particularly 
vivid example of a classical vortex which has many of the 
characteristics of a superfluid vortex has been realized in 
a recent experiment in which the 2n vortex phase winding 
was imprinted into a non-superfluid cloud of atoms and 
was observed to persist for extended times due to the 
cancellation of spins in the diffusion process [34[ ■ 

B. Off-diagonal long-range order 

Consider a system described by the usual one-body 
reduced density matrix 

p(r,r') = (* t (r)*(r')}, (2) 

where ^(r) is the second quantized bosonic field operator 
and the brackets denote quantum mechanical ensemble 
averaging. If p has a macroscopic eigenvalue No = O(N) 
where N is the number of particles in the system, it is 
said to possess Bose-Einstein condensation in the state 
determined by the corresponding eigenvector, also known 
as the macroscopic condensate wave function. 

The concept of long-range order is often encountered 
in the context of superfluidity. In a homogeneous Bose- 
Einstein condensed system, true long-range order exists 
in the sense that 

lim p(r, r') = const., (3) 

r— r' | — >oo 

where as in the normal state the off-diagonal correlations 
decay exponentially with the spatial separation 

p(r,r') ~ e -l r - r 'l/ So , (4) 
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where £o characterizes the length scale over which the 
correlations decay. Although the existence of Bosc- 
Einstein condensation does not a priori imply the sys- 
tem to be superfluid, the existence of off-diagonal order 
in the system can be considered as a prerequisite for su- 
perfluidity In the case of finite systems, such as trapped 
atomic gases, Eq. ([3]) generalises by understanding that 
the boundary of the system is mapped to infinity. 

C. Algebraic long-range order 



where K3 is the rate constant and rHot( r ) = (r)^(r)} 
is the total density, is an experimentally measurable 
quantity yielding overall information of the third order 
coherence properties of the system. While it has been 
used experimentally to infer the quasicondensation tran- 
sition point in helium thin films [231 ] , obtaining similar in- 
formation in trapped 2D gases is more problematic since 
the local gas density increases as the transition is crossed 
from normal to quasicondensed state and this partially 
compensates for the corresponding decrease in g( 3 \r). 



Two-dimensional systems lacking true long-range or- 
der may attain superfluidity through the Berezinskii- 
Kosterlitz-Thouless mechanism. In this case the corre- 
lations decay algebraically with distance 



p(r,r') 
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(5) 



and are characterised by the exponent, i](T) — n s fX^ B , 
where n s { is the 2D superfluid density, and AdB is the 
thermal de Broglie wavelength. For temperatures be- 
low the BKT transition temperature (Tbkt) the first 
order approximation to the critical exponent is i](T) = 
jT/Tbkt, corresponding to the universal jump at T = 
Tbkt in the superfluid density [7|]. 



D. Higher-order coherence 

Further insight into the state of the system is obtained 
by stud yin g the higher-order coherence properties of the 
system [35(. The second-order coherence function 



(2) (t)= (fr(r)frt( r )ft( r )ft( r )) 
" 1 (*t(r)*(r)) 2 



(6) 



yields information about the local coherence of the Bose 
field, ^, and in general for pth order coherence gc (r) = 
1 for n < p, while for the corresponding thermal state 
By decomposing the field into 'coherent' 



9t ( r ) 



— T7.I 



($) and 'incoherent' (ip) parts, * = $ + tp, Eq. © be- 
comes 
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where nth = fy'ip). In a similar fashion 
(|$| 6 + 9n th |$| 4 
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which essentially measures the probability of three- 
particle coincidences. The total three-body recombina- 
tion rate 



r = -K 3 I <7^(r)n tot (r) 3 d\ 



(9) 



E. Non-classical rotational inertia 

The moment of inertia, I(T), of a superfluid about 
a chosen axis is reduced from its classical value, I c \ = 
mN(r 2 ), due to the irrotational motion of superfluid mat- 
ter. The temperature dependent superfluid fraction 



N s f =1 _ J(T) 



Nt, 



Ic 



(10) 



where iV to t is the total particle number, is a quantity 
of special interest. A finite value of this macroscopic 
measurable may be used as an evidence of superfluid- 
ity. Microscopically, this information about superfluid- 
ity is encoded in the elementary excitation spectrum of 
the system. The collective scissors mode oscillation has 
been employed to prove that the occurrence of BEC in 
3D implies superfluidity [H, S3, H, !i| . Essentially, the 
scissors mode may be viewed as an oscillation of an ellip- 
soidal cloud of atoms about its semi-axis in the plane. In 
the collisionless regime, a gas in a normal state has two 
prominent undamped scissors mode eigenfrequencies 



LU± = \UJ X ± LU y \ 



(11) 



where uj x and u y are the planar trapping frequencies. In 
Eq. pip lo + corresponds to an irrotational quadrupole 
oscillation and a>_ is related to a classical rotational mo- 
tion. If superfluid component is present, it oscillates at 
an additional characteristic frequency 



(12) 



whose existence thus provides a clear sign for superflu- 
idity of the system in the collisionless regime. It is to 
be noted, however, that in the hydrodynamic limit ui- 
becomes over-damped and both the remaining thermal 
mode, oj+, and the superfluid scissors mode, w s f, attain 
the same value. In such situation, the damping rate of 
this mode may in principle be used to reveal superfluid re- 
sponse, although this may prove to be difficult to achieve 
in practice. 

The scissors mode excitations of the system are direct! 
related to the reduced moment of inertia 
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where Q(u>) is the Fourier transform of the time- 
dependent quadrupole moment Q(t) — J xy n to t( r )i) dr. 
Substitution of Eq. (fl~3| in Eq. (fTQ|) . yields a formula 
for the superfluid fraction in terms of the scissors mode 
excitations. We emphasize that the presence of a super- 
fluid scissors mode, Lo s f, implies non-classical rotational 
inertia, I < Z c j, finite superfluid fraction, N s f/N to t, and 
hence superfluidity. 



III. METHODS 

We employ the method of classical fields as detailed 
in Refs. SH, S3, El]- Essentially, this amounts to 
propagating the projected Gross-Pitaevskii equation 

iW t $> = -^V 2 $ + y cxt $ + 5 7>{|$| 2 $}, (14) 
Ira 

in time for the field, $(r, t), restricted in the subspace 
determined by an energy cut-off in the harmonic oscil- 
lator basis states. The projector, V, serves to constrain 
the evolution of the field within the subspace of highly 
occupied states. Each simulation corresponds to the evo- 
lution of a single trajectory through the phase space and 
therefore, in order to construct thermodynamic observ- 
ables, one should ensemble average over many different 
but equivalent trajectories. However, when consider- 
ing an equilibrium quantities, we may assume the sys- 
tem to be ergodic and replace such ensemble averages 
by time-averages over the instantaneous field configura- 
tions taken from a single trajectory. The field, <E>(r, t), 
is normalized to the number of particles, N c i, described 
by the restricted basis. The total number of particles, 
Not = N c i + Nth, is obtained by using the semiclassi- 
cal Hartree-Fock approximation for the Nth above cut-off 
particles as in Refs. [lj| H|| and as described below. The 
in-plane phase function, ip(x,y,t), of the complex field, 
<&(r, i), allows for an explicit detection of the locations of 
vortices and antivortices. The classical field is completely 
described by the conserved total energy, E c \, an energy 
cut-off for the restricted basis, E cut , the dimensionless 
nonlinearity constant, C = gN c \/hu x a x , and the har- 
monic trap frequencies, lo x ,u> v and u) z . Here the spatial 
length scale is a x = \J h/2muj x . From these simulations 
we can also compute the equilibrium temperature, T, and 
chemical potential, fi, as described in Refs. [iH |46|. 



of the lowest axial oscillator states, n z = {0,1,2...}, 
contribute significantly to the total number of particles. 
However, the temperature is too low for the equipartition 
theorem to apply and therefore these axial levels need to 
be treated discretely in the semiclassical integral. The 
Hartree-Fock energy for this system is given by 

E^(K, x, y)=K + V CKt {x, y, 0) + (n z + l/2)faj z 

+ 2^ff™,,rXh(^)+ 2 5".,0<l$(^y)l) 2 , (15) 
m 

with K the kinetic energy. The thermal densities are 
given by 
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e {E^{K,x,y)-n)/k B T _ l 



dK, 



where 



(16) 



K min = max{0, E cut - V cxt (x, y, 0) - fru z /2}, (17) 

and n is the chemical potential. The interaction term in 
Eq. (Tl5|) contains the multilevel coupling constant 



(18) 



9 I \<t>n{z)\ \<l>m{z)\ dz, 



which accounts for the interactions between the particles 
in different axial levels. We have denoted y, z) — 
<$>(x,y)(j)o(z) and (f> n (z) are the normalized harmonic os- 
cillator eigenstates. Finally, the number of above cut-off 
particles 



Nth = J n th( x >y) dxd vi 



(19) 



is obtained by integrating over the 2D densities and sum- 
ming the contributions from different axial energy levels. 
To obtain data points for a fixed number of particles AT tot 
for a range of temperatures, we estimate the cutoff en- 
ergy -E cu t, total number of particles N c \ and energy E c \ 
of the classical re gion of the system according to the pre- 
scription of Ref. [43] . We then simulate this within the 
PGPE and calculate the temperature and total number 
of thermal particles using the procedure described above. 
We use this knowledge to make any necessary adjustment 
to the initial guesses to end up with a target iVtot- 



Collective excitations 



A. Semiclassical approximation 

An inherent feature in our numerical method requires 
that the total particle number, Atot, must be computed 
a posteriori. This is done within the self-consistent 
Hartree-Fock approximation by computing the particles 
not included in the simulated field, $, from the semi- 
classical density. The form of the semiclassical integral 
reflects the quasi-2D nature of the trap. For the temper- 
atures considered here, ksT ~ hw z and therefore several 



In principle we could construct an approximation to 
the full Green's function from our time-dependent clas- 
sical field simulations allowing us to extract the col- 
lective excitation frequencies and their damping rates 
for the system. It would be, however, a formidable 
task. Furthermore, an identification of near-degenerate 
modes would become cumbersome. Instead, we concen- 
trate on a specific class of excitations — the so-called scis- 
sors modes — which may be selectively excited and con- 
sequently their oscillation frequency can be individually 
measured from our dynamical simulations. 
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In order to accurately compute the scissors mode col- 
lective oscillation frequencies, an ensemble averaging over 
many equivalent trajectories is required. Therefore we 
first prepare a large set of initial field configurations by 
time-sampling a single equilibrium simulation. These in- 
stantaneous field configurations are then rotated 11 de- 
grees in the x — y -plane with respect to the semi-axis of 
the anisotropic trapping potential. To facilitate lossless 
rotation of the state, we must first project the classical 
field into a larger eigenbasis in order to account for the in- 
crease in the energy after rotation due to the anisotropic 
trapping geometry. Subsequently the seed states thus 
prepared are propagated in time and the quadrupolc mo- 
ment, Qi{t), is measured. The above procedure is re- 
peated for S = 150 — 200 microstates for each temper- 
ature point. The scissors mode frequencies are then ob- 
tained from the Fourier transformation of the ensemble 
averaged quadrupole signal Q(t) — J2i =1 Qi(t)/ S , from 
which also the damping rates of those modes can be es- 
timated. 



C. Computational parameters 

In order to allow practical comparison with future ex- 
periments, we choose experimentally realistic system pa- 
rameters which allow us to study the scissors-mode col- 
lective oscillation frequencies. The strength of particle 
interactions is determined by a constant, g = 4Trh 2 a/m, 
and the harmonic confining potential, V cx t = m(uj 2 x 2 + 
WyU 2 + w5rz 2 )/2, is characterized by the Cartesian fre- 
quencies {uj x , uj y ,u! z } — 2tt x {10, 20, 4000} Hz. The trap 
is chosen to be anisotropic in the x — y -plane in order to 
lift the degeneracy between the quadrupole and scissors 
modes. However, the planar anisotropy is kept moder- 
ate, in contrast to Ref. 29], in order to separate the fre- 
quency of the superfluid scissors-mode from both of the 
two classical modes. We consider AT tot = 10 5 87 Rb atoms 
interacting with the s-wave scattering length a = 5.3 nm. 



D. Radial averaging 

Several results we present are generated by performing 
an average along elliptical trajectories in the 2D plane 
about the trap centre. This is done in order to utilize 
the full information contained in the simulated fields. 
Since our system does not possess cylindrical symmetry 
(in order to facilitate computation of scissors modes), we 
perform this averaging by considering elliptical shells of 
constant (2D) trap potential energy and average over all 
spatial points falling on such strips. In what follows, ro, 
denotes the distance from origin to such ellipse along the 
weakly trapped x-axis of the trap. 



IV. RESULTS 

In this section we present our numerical results and 
analysis for the superfluid indicators described above. 
Taking the compendium of results described below, we 
claim that the quasicondensate studied here is superfluid. 
All measurables extracted from our simulations point to 
a cross-over temperature T s f at which our system attains 
superfluidity. 

A. Fluctuations and vortices 

In order to lay down the qualitative features of the 
system we have plotted instantaneous 2D classical field 
densities, X 2 lB \^>(x,y)\ 2 , in Fig.QJa) and (b), and the cor- 
responding phases ip(x,y), Fig. QJc) and (d). The tem- 
peratures are T = 114 nK for Fig. [IJa) and (c) and 
T = 151 nK for Fig. QJb) and (d). At low tempera- 
tures the density and phase are relatively uniform, while 
at high temperatures both exhibit strong fluctuation and 
and vortices and antivortices are nucleated. This obser- 
vation is in striking contrast to the usual situation in 3D 
and highlights the main qualitative difference between 
2D and 3D systems. 

To further quantify the emergence of vortices and an- 
tivortices due to the phase fluctuations, we have mea- 
sured at each temperature point the probablility, P v (r ) 1 
of finding a vortex or antivortex at radius tq. This is 
done by locating all phase singularities in an instan- 
taneous classical field configuration and averaging over 
1000 different microstates. The classical field area is 
then divided into ellipsoidal strips of equal width and 
the vortex occupation probability is obtained by count- 
ing the number of phase singularities detected within 
each strip divided by the number of microstates sam- 
pled (this is the radial averaging discussed earlier) . Thus 
obtained probabilities for a set of temperatures, {Tj} = 
{168, 158, 155, 148, 143, 138, 135, 131, 126, 120, 115} nK, 
are plotted as functions of radial distance, ro, in Fig. [5] 
The bullets denote the coherence length, measured as 1 /e 
radius of g(0, r ) (defined in Sec. C below) and the red 
curve is for T s f = 155nK. There is a sudden jump in 
the vortex occupation probability at the cross-over. In 
the superfluid phase there is a vortex-free region at small 
radii. Vortex pairs are observed in a narrow band near 
the edge of the coherent region of the system, i.e. the 
spatial region that has a flat phase in FigQJc). Thus the 
system can be divided in three concentric regions: central 
coherent and vortex free BEC-like region, coherent BKT- 
like region where vortices are bound, and an incoherent 
thermal outmost region where vortices are free. 

B. Condensate fraction 

We compute the one-body density matrix, Eq. (TTJ, for 
the classical field by assuming ergodicity, which allows 
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FIG. 1: (Color online) Density (a),(b) and phase (c),(d) of the classical field at two different temperatures: (a),(c) T — 114 nK, 
and (b),(d) T = 151 nK. Vortices and antivortices are denoted by + and — signs respectively. 




20 

r (/j,m) 

FIG. 2: (Color online) Radial vortex occupation probability 
density at different temperatures (from left to right) {Ti} — 
{168,158,155,148,143,138,135,131,126,120,115} nK. The 
red line indicates the superfluid crossover at T a f = 155 nK. 



with the pure 2D ideal gas relation 



N \T ) ' 



(20) 



where To is the critical temperature calculated for 
a quasi-2D ideal-gas in our trap geometry containing 
Ntot = 1.2 x 10 5 particles. The vertical line is the temper- 
ature T s f = 155 nK, below which this system is superfluid 
and the markers are the simulation data. Unlike in the 
recent experiment of Kriiger et. al [30 ] we find the inter- 
actions to cause only a minor shift in the critical temper- 
ature from the non- interacting boson result. In terms of 
the definition based on the eigenvalues of the density ma- 
trix, the system may be claimed to show Bose-Einstein 
condensation at all temperatures below T s {, Nevertheless 
it turns out that the system is best described in terms 
of phase fluctuating quasicondensate apart from the very 
lowest temperatures. 



Coherence 



us to replace the ensemble average by a time average. 
The number of condensed particles, N , is obtained by 
computing the largest eigenvalue of the density matrix. 
Figure [3] displays the condensate fraction as a function of 
temperature. The curve is plotted to provide comparison 



Coherence is an essential feature of superfluidity. The 
density matrix, Eq. (Q}, provides an useful probe for the 
global coherence between two spatially separated points 
in the system. In particular it conveys the knowledge of 
the correlation length and the information on the possible 
presence of long-range order. Figure H shows the two- 
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T (nK) 

FIG. 3: (Color online) Condensate fraction as determined by 
diagonalization of the one-body density matrix as function of 
temperature. The markers are the data, the black curve is 
the 2D ideal gas result and the red vertical line indicates the 
superfluid crossover temperature T s f. 



point function, 

a(0r)= ($*(0)$M) 

' V(l'i > (0)| 2 +«th(0))(|$(r )|2+n th (ro)) 

for different temperatures. For low temperatures, T < 
T s { , and small radii we witness power-law decay of g(0, ro) 
in accordance with Eq. ([5]), where as for temperatures, 
T > T s f, and/or near the edge of the coherent region, 
exponential decay is observed. The qualitative behavior 
changes at the cross-over temperature, T s f, denoted by 
the red line. 

The second-order coherence function, Eq. mea- 
sures local coherence in the gas. Particularly, for a purely 

(2) 

thermal sample g\ (r) = 2 and for completely coherent 

state g^c\v) — 1. In our inhomogeneous system g^ 2 \r) 
interpolates between these two values as shown in Fig. [5] 
where g^(ro) is displayed for different temperatures as 
functions of the radial distance from the trap centre. At 
the lowest temperatures, <?' 2 '(ro), shows a flat part in 
accordance with the presence of a nearly phase coherent 
BEC. At large enough radii and for high enough temper- 
atures g^ (r ) approaches its thermal value. The plateau 
disappears near the temperature where the macroscopic 
wave function vanishes. The thicker red line indicates 
g {2) (r ) for a result with T ~ T sf . 

In Fig. [6l we have plotted the coherence length as func- 
tion of temperature. The values are obtained by measur- 
ing the 1/e width of the two-point correlation function. 
We have also plotted a function, £-(T) = \t\ p , where, 
t oc (T sf - T)/T sU (solid curve) for p = 0.25 and T < T sf . 
The value for the exponent, p, may be crudely explained 

1 /4 

in terms of the Thomas- Fermi radius, Rtf oc N q of 
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FIG. 4: (Color online) Two-point correlation functions, 
3(0, ro) as function of spatial distance, ro, from the trap cen- 
tre for a range of temperatures, {Ti}. The thicker red line is 
for the estimated superfluid cross-over temperature. 




1 10 20 30 40 

r (/im) 

FIG. 5: (Color online) Second-order coherence functions 
9 ( r o) as function of spatial distance, ro from the trap cen- 
tre for a range of temperatures, {Ti}. The thicker line is for 
the estimated superfluid cross-over temperature. 



an isotropic 2D condensate, since N /N tot varies linearly 
with the temperature in the vicinity of the cross-over 
point, T sf . 



D. Scissors modes 

We have computed the scissors mode frequencies ac- 
cording to the description in the methods section. The 
obtained oscillation frequencies and their relative inten- 
sities are displayed in Fig. [71 The blue bullets are the 
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FIG. 6: (Color online) Correlation length as a function of 
temperature. The bullets are the numerical data, the vertical 
red line denotes the superfiuid cross-over temperature and the 
black curve is the power-law function described in the text. 



FIG. 7: (Color online) Scissors mode frequencies as a function 
of temperature. The horizontal dashed lines are the analytical 
predictions of Eqs. (|lip and (|12|l and the solid vertical line 
is our estimate of the superfiuid transition temperature. The 
horizontal axis refers to the equilibrium temperature prior to 
the rotation of the state. The normalised Fourier spectrum 
at each temperature is indicated in gray in the background. 



mean oscillation frequencies obtained by fitting double 
Gaussian functions to the normalised Fourier spectrum 
at each temperature, which is indicated in gray in the 
background. The horizontal dashed red lines are the an- 
alytical predictions Eqs. (fTTj) and (fT2")l . Above the cross- 
over temperature, T s f, indicated by the red vertical line, 
we obtain signal for two thermal modes whose frequencies 
are found to agree with the predictions of Eq. (fTTj) . At 
temperatures well below the cross-over, only one scissors 
mode persists with a frequency corresponding to that of a 
superfiuid, see Eq. (fl2|) . This is the key feature verifying 
the quasicondensate to be superfiuid. 

It is interesting to notice that we only observe two 
different scissors modes at all temperatures. The up- 
per 'irrotational' mode simply experiences frequency shift 
across the cross-over associated with the change in su- 
perfiuid density of the system. This is to be contrasted 
with 3D systems where in general three different scis- 
sors modes exist and the superfiuid scissors mode experi- 
ences downward (as opposed to the upward shift seen in 
Fig. [7]) frequency-shift on increasing temperature across 
the cross-over [37|. 

We have also computed the scissors modes for different 
systems of varying particle interaction strengths which 
could be experimentally realized using Feshbach reso- 
nances. We have verified that in the strongly interacting 
systems only a single scissors mode survives at all tem- 
peratures, making it difficult to distinguish the superfiuid 
and thermal response of the system from one another. 



V. DISCUSSION 

In conclusion, we have studied the problem of 
superfluid-normal cross-over in a real, trapped quasi- 
2D Bose gas. In such systems the formation of pure 
Bose-Einstein condensation is challenged by the long- 
wavelength phase fluctuations and this fact has made 
the characterization of such systems difficult both the- 
oretically and experimentally. By performing classical 
field simulations for these systems, we have shown that 
such quasicondensates are superfiuid below the cross- 
over temperature T s {. This conclusion is based on ob- 
servations of the global coherence properties and scis- 
sors mode excitations of the system, which constitute 
the two major results of this article proving superfluid- 
ity of quasicondensates. In particular, the emergence of 
condensate scissors-mode below the cross-over temper- 
ature provides an unequivocal microscopic evidence of 
non-classical rotational inertia and thus the superfluidity 
of the system. We have not found signs of fragmentation 
in terms of the eigenvalues of the density matrix below 
the cross-over temperature. A similar conclusion is ob- 
tained from the fact that only a single condensate peak 
is observed in momentum space. This indicates that de- 
spite of the prevailing phase fluctuations, the superfiuid 
state of these systems resembles more closely that of a 
single-mode trapped Bose-Einstein condensate than the 
bulk Berczinskii-Kosterlitz-Thoules superfiuid phase. 

Quasi- 2D quantum gases are currently under active ex- 
perimental investigation. While experiments have veri- 
fied aspects of first-order coherence in this system consis- 
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tent with the Berezinskii-Kosterlitz-Thouless transition, 
there has yet to be any direct evidence of superfluidity. 
In this paper we have reported a series of tests that can 
be realized in an experiment and should cast light on the 
formation of superfluidity in these systems. 
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